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20.  Abstract 

thermohydrodynamic  structure  of  each  vortex  Is  being  developed,  and  then  the 
conditions  for  transition  from  the  moderately  intense  to  the  very  Intense 
vortex  are  to  be  sought.  The  practical  motivation  Is  to  make  progress  toward 
the  highly  desirable,  but  very  forml dab! e> task  of  being  able  to  anticipate 
which  tropical  storms  or  minimal  hurricanes  will  evolve  to  supertyphoons. 

Since  the  one-cell  storm  Is  believed  to  be  largely  described  already  for 
present  purposes,  description  of  the  two-cell  vortex  Is  the  current  challenge 
In  particular,  the^bropertles  and  location  of  the  eye  wall  are  sought,  since 
the  potential -vortex  and  surface-inflow  subdivisions  of  the  structure  of  a 
very  intense  vortex  are  In  hand. 

The  modeling  proceeds  from  basic  thermohydrodynamic  principles.  Non- 
essential  geometric  detail,  as  well  as  association  of  conclusions  with  the 
details  of  particular  parameterlzatlons ,  Is  being  avoided  as  much  as  possible 
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PREFACE 
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investigator,  and  Phillip  Feldman,  numerical  analyst,  of  TRW  Defense  and 
Space  Systems  Group,  and  George  Carrier  of  Harvard  University,  Cambridge, 

MA,  consultant  to  TRW.  Also  contributing  was  Paul  Dergarabedian  of  the 
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Martin  and  James  Hughes  of  the  Atmospheric  Sciences  Program  of  the  Office  of 
Naval  Research  for  their  encouragement  and  their  patience. 
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SUMMARY 


Analytic  modeling  of  well -organi zed  rotating  convective  storms  is 
being  taken  to  elucidate  the  evolution  from  a  moderately  intense  one-cell 
vortex,  characterized  by  low-level  pressure  deficits  on  the  order  of  one 
percent  of  atmospheric  pressure,  to  a  very  intense  two-cell  vortex,  char¬ 
acterized  by  low-level  pressure  deficits  on  the  order  of  ten  percent  of 
atmospheric  pressure.  The  physical  distinction  between  the  two  stages  is 
the  insertion  of  a  dry,  compressionally  heated,  nonrotating,  central  down- 
draft  of  originally  tropopause-level  air  in  the  more  severe  case.  The  quasi¬ 
steady  mature  description  of  the  thermohydrodynamic  structure  of  each  vortex 
is  being  developed,  and  then  the  conditions  for  transition  from  the  moderately 
intense  to  the  very  intense  vortex  are  to  be  sought.  The  practical  motivation 
is  to  make  progress  toward  the  highly  desirable,  but  very  formidable 
task  of  being  able  to  anticipate  which  tropical  storms  or  minimal  hurricanes 
will  evolve  to  supertyphoons. 

Since  the  one-cell  storm  is  believed  to  be  largely  described  already 
for  present  purposes,  description  of  the  two-cell  vortex  is  the  current 
challenge.  In  particular,  the  properties  and  location  of  the  eye  wall  are 
sought,  since  the  potential -vortex  and  surface-inflow  subdivisions  of  the 
structure  of  a  very  intense  vortex  are  in  hand. 

The  modeling  proceeds  from  basic  thermohydrodynamic  principles. 
Nonessential  geometric  detail,  as  well  as  association  of  conclusions  with 
the  details  of  particular  parameter izations ,  is  being  avoided  as  much  as 
possible. 

I.  INTRODUCTION 

(a)  STATEMENT  OF  THE  OBJECTIVE 

The  principal  contribution  still  to  be  made  in  forecasting  intensity 
of  tropical  storms  is  to  identify  the  observable(s)  that  permit  one  to 
anticipate  (as  much  in  advance  as  possible)  which  incipient  storms  will 
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become  hurricanes  and  which  will  not. 

Very  roughly  about  half  of  all  tropical  storms,  one-cell  vortices 
extending  vertically  from  sea  level  to  torpopause  and  radially  many  hundreds 
of  miles,  transform  into  two-cell  vortices  (hurricanes)  (Fendell  1974). 

The  transition  involves  insertion  of  a  column  of  relatively  dry,  clear, 
nonrotating,  originally  tropopause-1 evel  ,  compressionally  heated,  slowly 
recirculated  air  at  the  center  of  the  vortex.  Accordingly,  the  region  of 
intensely  swirling,  very  cloudy,  torrentially  raining  updraft  becomes 
displaced  to  an  annulus.  An  "eye"  has  then  been  inserted  within  an  "eyewall". 

A  formidable  vortex,  with  peak  swirl  of  perhaps  not  100  mph  and  horizontal 
pressure  deficit  at  sea  level  of  a  few  tens  of  millibars  (the  tropical  I 

storm,  as  defined  here),  has  transformed  into  a  terrifying  vortex  of  possibly 
over  200  mph  peak  swirl  and  a  pressure  deficit  of  100  mb  or  more  (the 
hurricane,  as  defined  here).  The  question  arises,  what  telltale  observable 
accessible  to  aircraft,  satellites,  etc.,  would  permit  one  to  anticipate 
the  transition  from  tropical  storm  to  hurricane?  Of  course,  the  intensity 
of  a  tropical  cyclone  is  not  a  simple  monotonic  increase  with  time  to  a  peak, 
following  by  a  simple  monotonic  decay;  thus,  there  may  be  several  transitions 
back  and  forth  between  one-  and  two-cell  structure,  suggestive  of  partial 
insertion,  then  removal,  of  an  “eye"  (Carrier  1971a). 

The  initial  goal  of  this  project  is  the  formidable  one  of  delineating 
the  structure  of  a  well -developed  two-cell  vortex,  to  complement  the 
relatively  wel 1 -understood  structure  of  a  one-cell  vortex.  Before  subleties 
of  the  transition  can  be  detected,  the  "stable  equilibrium  states"  should 
be  characterized  in  sufficient  quantitative  detail. 

Normally  quantifying  any  details  about  structure  in  tropical  cyclones 

#  Gray  (1979,  p.  51)  comments:  "Tropical  cyclone  forecasters  generally 
agree  that  skill  at  operational  forecasting  of  tropical  intensity  change 
is  nearly  zero."  He  notes  that  trying  to  relate  intensity  to  cloudiness 
configuration  (as  observed  by  passive  satellites)  is  subject  to  signifi¬ 
cantly  erroneous  results,  and  furnishes  inadequate  detail.  It  may  be  noted 
that  the  gustiness  of  winds  is  not  discussed,  since  treatment  of  the 
variance  seems  premature  when  prediction  of  the  mean  is  still  not  possible. 

Gustiness  factors  for  tropical  cyclones  (ratio  of  peak  transient  wind  to 
wind  sustained  over  a  scale  of  one  minute  to  five  minutes,  depending  on 
the  correlator!  range  from  1.0  to  1.8  in  the  literature  (Brand  et  al  1979). 
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is  highly  uncertain  because  proper  formulation  of  turbulent  diffusion, 
cummulus  convection,  and  radiative  transfer  is  unknown.  These  three 
important  small-scale  phenomena  may  well  not  have  satisfactory  representation 
in  terms  of  macroscale  variables,  so  even  the  ve-y  attempt  at  such 
parameterization  is  doomed  at  the  outset.  Yet  such  parameterization  is 
the  essence  of  all  contemporary  tropical  cyclone  modeling,  and  any  new 
undertaking  could  just  augment  the  plethora  of  almost  totally  uncorroborated 
speculation.  Why  is  there  cause  for  optimism? 

(b)  AN  APPROACH  TO  TWO-CELL  STRUCTURE  FOR  AN  INTENSE  ATMOSPHERE  VORTEX 

The  answer  concerns  a  little  ingenuity,  a  little  insight,  and  a  lot 
of  luck. 

The  ingenuity  lies  in  the  fact  that,  whereas  almost  all  contemporary 
hurricane  modeling  entails  direct  numerical  assault  on  all-inclusive 
mathematical  formulation,  this  project  utilizes  an  alternative,  approximate 
analytic  treatment  as  much  as  possible.  While  other  liabilities  may  be 
cited  (e.g.,  finite-differencing  introduces  spurious  effects  not  present 
in  the  physical  formulation),  there  is  one  particular  difficulty  with 
direct  numerical  computation  of  special  relevance  here:  a  physical  effect 
important  anywhere  in  the  flow  field  must  be  retained  everywhere  in  the 
flow  field.  Thus  rather  arbitrary  models  for  turbulent  diffusion,  cumulus 
convection,  and  radiative  transfer  are  present,  possibly  to  vitate  results 
throughout  the  flow  field.  In  the  alternative  approximate  aralytic  approach 
examined  here  (Carrier  1970;  Carrier,  Hammond  and  George  1971),  the  hurricane 
is  subdivided  logically  into  those  parts  in  which  different  physical 
phenomena  dominate;  when  a  particular  subdivision  is  treated,  only  that 
subset  of  the  full  equations  that  retains  the  locally  dominant  physical 
processes  is  solved.  Such  a  subset  is  usually  more  tractable  to  solve,  and 
parametric  variation  easier  to  examine.  A  composite  global  solution  is 
synthesized  from  solutions  to  the  subdomains  by  demanding  appropriate 
continuity  of  dependent  variables  and  their  fluxes  at  interfaces. 

The  little  insight  referred  to  above  is  the  deduction  that,  under  the 
high-speed  portion  (only)  of  a  severe  rotating  storm,  the  nonlinear  "surface" 
Inflow  layer  is  inviscidly  controlled;  only  in  a  very  small  sublayer 
(immediately  continguous  to  the  ground)  of  the  inflow  layer  does  diffusion 
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enforce  the  no-slip  boundary  condition  (Carrier  1971b;  Burggraf,  Stewartson 
and  Belcher  1971;  Carrier  and  Fendell  1978).  The  diffusive  sublayer 
becomes  thinner,  and  the  outer  inflow  layer  thicker,  as  one  moves  in  toward 
the  center  of  the  vortex.  Why  this  is  so  is  subtle,  but  it  is  familiar 
in  nonrotating  contexts  to  aerodynamicists  in  the  form  of  the  following  bit 
of  well-known  empiricism  (e.g..  Launder  1964):  a  near-wall  diffusive  layer 
becomes  thinner  in  the  direction  of  a  favorable  (i.e.,  accelerating) 
pressure  gradient.  At  sea  level,  the  pressure  is  high  at  the  lateral 
edge  of  an  intense  vortex  and  the  pressure  is  greatly  reduced  at  the  center, 
such  that  a  strong  positive  pressure  gradients  exists;  thus,  the  low-level 
swirling  influx  (that  erupts  up  the  eyewall )  is  predominantly  inviscidly 
controlled  as  it  proceeds  toward  the  hurricane  center. 

Now,  radiative  cooling  is  important  only  in  the  outer  and  less  cloudy 
portions  of  the  storm  (Fendell  1974),  and  the  simple  well-defined  model  of 
moist-adiabatic  ascent*  suffices  to  describe  the  locus  of  thermodynamic 
states  in  the  eyewall.  In  the  previous  paragraph  it  is  pointed  out  that 
the  low-level  flow  entering  the  eyewall  is  nondi ffusive .  This  is  the  crucial 
ingredient  of  good  fortune  that  joins  ingenuity  and  insight.  Two  of  the 
three  small-scale  processes  that  are  the  bane  of  the  numerical  modeler 
(specifically,  turbulent  diffusion  and  radiative  transfer)  do  not  enter 
significantly  in  the  portion  of  the  hurricane  flow  field  critical  to 
quantitative  analysis  of  two-cell  structure,  and  the  third  process  (cumulus 
convection)  can  be  formulated  locally  in  a  well -accepted ,  well-defined  form. 

The  upshot  is  the  following  good  fortune:  one  of  the  few  questions 
about  hurricane  structure  that  can  be  definitively  formulated  by  contemporary 
modeling  happens  to  be  one  of  the  most  important.  Specifically,  that 
question  is,  what  is  the  balance  of  forces  and  thermodynamic  states  in 
the  central  region  of  a  tropical  storm  that  permits  two-cell  structure 
to  arise  as  an  alternative  to  one-cell  structure? 

*  Water-vapor-containing  air  rises,  expands  and  cools  dry-adiabatical ly 
unci  1  saturation.  Thenceforth,  it  continues  to  ascend  such  that  saturation 
of  water  vapor  is  maintained  at  the  local  thermodynamic  state,  the  condensed 
excess  water  substance  falling  out  so  that  it  is  unavailable  for  re-evapora¬ 
tive  cooling.  (Hence  the  process  cannot  be  reversed.)  The  originally  low- 
level,  convectively  unstable  fluid  undergoing  this  ascent  is  taken  to  rise 
too  quickly  to  entrain  any  ambient  air,  and  hence  remains  unmixed  with  cooler 
or  drier  air.  In  fact,  of  course,  some  mixing  with  nearby  ambient  air  occurs 
and  some  re-evaporation  into  this  unsaturated  mix  also  occurs,  but  the  no- 
mixing  case  is  realistic  enough  to  provide  valuable  insight. 


(c)  THE  BOUNDARY-VALUE  PROBLEM 


For  axi symmetric  inviscid  swirling  flow  in  which  a  radial  influx 
becomes  an  axial  upflux,  the  boundary-value  problem  is  well-defined. 
Especially  if  one  confines  attention  at  the  outset  to  low  altitudes  such 
that  incompressibil ity  suffices,  the  formulation  is  available  in  texts 
(e.g.,  Batchelor  1967).  Angular  momentum  is  constant  on  streamlines, 
continuity  is  enforced,  advective  and  centripetal  accelerations  balance 
the  radial  pressure  gradient,  and  the  total  pressure  head  is  constant 
across  streamlines  (Carrier,  Dergarabedian  and  Fendell  1969).  The  annulus 
is  bounded  on  one  side  by  a  free  streamsurface  that  demarcates  the  eyewall 
interface  with  the  eye,  and  on  the  other  side  by  a  free  streamsurface 
that  demarcates  the  eyewall  interface  with  the  outer  vortex.  The  position 
of  both  these  bounding  streamsurfaces  must  be  found  in  the  course  of 
solution,  under  enforcement  of  statements  concerning  the  continuity  of 
pressure.  One  result  sought  is  the  displacement  of  the  eyewall  from  the 
axis  of  symmetry,  which  is  also  the  axis  of  rotation.  It  is  anticipated 
that  this  displacement  ultimately  increases  with  height. #  The  mathematically 
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The  increasing  displacement  of  the  eyewall  from  the  axis  with  altitude 
seems  likely  under  conservation  of  angular  momentum,  since  density 
decreases  with  altitude.  Conservation  of  mass  also  suggests  that  the 
cross-sectional  area  of  the  eyewall  annulus  increases  with  altitude.  The 
fact  that  the  eyewall  tilts  radially  outward  with  increasing  height  is 
suggested  further  by  certain  satellite  photographs  of  'hurricanes 
(Fendell  1974).  This  outward-sloping  behavior  may  seem  just  another 
passing  observation,  but  actually  it  bears  directly  on  perhaps  the  most 
controversial  question  dividing  hurricane  models  during  the  last  three 
decades.  That  question  concerns  the  amount  of  augmented  sensible  and 
latent  heat  transfer  (above  ambient  tropical  autumnal  transfer  rates) 
from  ocean  to  atmosphere  that  is  necessary  near  the  center  of  a  hurricane 
to  sustain  the  hurricane  (Carrier,  Hammond  and  George  1971).  If  the 
eyewall  were  vertical,  augmentation  of  enthalpy  transfer  from  sea  to 
air  would  have  to  be  an  order  of  magnitude  greater  than  conventional 
rates  (Malkus  and  Riehl  I960);  pressure  deficits  achievable  under 
compressional  heating  in  the  eye  could  not  serve  as  an  explanation  of 
intense  swirl  speeds  observed  in  the  eyewall,  and  instead  oceanic  heat 
transfer  would  have  to  explain  the  density  reduction  that  hydrostatically 
integrates  to  low  pressure  at  the  base  of  the  eyewall.  Such  an  increase 
in  oceanic  transfer  seems  mechanistically  unlikely.  lf_  the  eyewall 
slopes  radially  outward,  no  such  augmentation  in  enthalpy  transfer  is 
requisite,  because  low  pressure  at  the  base  of  the  eyewall  could  be 
explained  by  density  reduction  through  compressional  heating  in  the 
upper  troposphere. 
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elliptic  character  of  the  boundary-value  problem  reflects  the  physical 
fact  that  streamwise  and  transverse  derivatives,  and  streamwise  and  trans¬ 
verse  velocity  components,  are  of  comparable  magnitude;  there  is  no  single 
preferred  direction  of  change.  Indeed,  frankly,  a  great  deal  of  time  has 
been  unsuccessfully  expended  by  the  investigators  seeking  a  valid  mathematical 
transformation  that  would  convert  the  elliptic  problem  into  a  parabolic 
one;  an  axisymmetric  parabolic  problem  is  more  readily  solved  precisely 
because  it  possesses  a  single  dominant  direction  of  change.  Such  a 
transformation  seemed  attained,  but  it  ultimately  proved  a  chimera,  and 
now  the  investigators  must  come  to  grips  with  the  nonlinear  elliptic 
probl em . 

Some  preliminary  insight  is  attainable  at  the  sacrifice  of  detail  by 
using  an  "integral  method"  (Finlayson  1972).  Roughly,  one  "forces"  a 
parabol ic-1  i ke character  on  the  problem  by  taking  transverse  averages  and 
then  seeking  streamwise  variance  of  these  averages  from  an  appropriately 
simplified  boundary-value  problem.  Such  gross  resolution  can  suffice  as 
preliminary  guidance  only,  and  ultimately  greater  resolution  is  required 
for  present  purposes. 

(d)  ON  SIMILARITIES  OF  SWIRL-STA8 II I  ZED-COMBUSTOR  FLOWS  AND  HURRICANE 

DYNAMICS 

The  meteorologically  motivated  swirling-flow  studies  being  discussed 
should  be  of  value  in  aerodynamically  stabilized  continuous  combustors. 

Such  combustors  are  used  in  face-fired  furnaces  (industrial  boilers,  whether 
fueled  by  oil  or  coal)  and  in  gas  turbines.  In  fact,  the  dual-purpose 
pertinence  of  the  proposed  research  is  such  that  the  matter  is  now  developed 
at  some  length. 

For  fuel  efficiency,  low  pollution,  and  flame  stabilization,  combustors 
with  recircul atory  zones  are  quite  common  (Beer  and  Chigier  1972).  For 
premixed  fuel  and  air,  such  a  reversed-flow  domain  lets  one  flow  the  fresh 
charge  at  speeds  in  excess  of  the  flame  speed,  without  having  blow-off. 

For  diffusion  flames  (unpremixed  fuel  and  air),  one  can  retain  the  reactants 
longer  in  the  combustion  zone  if  either  reduced  volatility  or  slow  gas-phase 
reaction  rates  (relative  to  flow  rates)  poses  a  problem.  One  may  view  the 
mechanism  by  which  the  recirculation  zone  succeeds  in  bringing  continuously 
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introduced  fresh  charge  to  its  ignition  condition  in  several  ways.  The 
recirculatory  "bubble"  retains  heat  and  hot  radicals  to  raise  the  charge 
temperautre  and  to  instigate  chain-branching  chemical  mechanisms.  Also, 
the  stagnation-point-generating  recirculatory  "bubble"  reduces  flow  speed 
relative  to  reaction  rate,  with  the  result  that  less  incompl etely-burned 
reactant  escapes;  such  "lost"  reactant  is  not  only  economically  wasteful 
and  polluting,  but  energetically  unfavorable,  because  some  sensible  heat 
is  convected  away  without  exploitation  of  the  full  chemical  exothermic i ty 
possible. 

Recirculation  may  be  achieved  in  the  wakes  of  bluff  bodies  introduced 
into  the  flow,  but  these  entail  large  pressure  losses  in  high-speed  flow, 
and  in  general  introduce  material s-survival  and  coke-deposition  problems. 

In  recent  years,  abrupt  area  expansion  of  the  flow,  secondary  air  jets, 
and  other  adverse-pressure-gradient-generating  designs  have  been  implemented 
to  achieve  recirulatory  flow  aerodynamical ly,  i.e.  without  introduction 
of  physical  obstacles  into  the  flow  (Syred,  Chigier  and  Beer  1971).  However, 
of  prinicpal  interest  here,  is  the  use  of  swirl ing  flow  at  the  inlet,  either 
of  both  the  fuel  and  air  stream,  or  at  least  of  the  air  stream  alone,  to 
aid  achievement  of  recirculatory  flow  in  the  reaction  zone  of  the  combustor. 

In  general  in  swirling-stabilized  combustors,  the  reactants  enter  the 
combustor  with  axial  (i.e.,  streamwise)  momentum  and  with  angular  momentum; 
the  ratio  of  the  latter  to  the  former  is  defined  as  the  swirl  ratio  and 
traditionally  denoted  S  (Beer  and  Chigier  1972).  [Some  authors  (e.g.. 

Hall  1972)  define  an  equivalent  angle  of  swirl,  $,  by  the  definition  tan;  s 
(v/w),  where  v  denotes  the  swirl  component  of  the  velocity,  and  w  denotes 
the  axial  component.]  Typically,  for  S  <  0.6,  the  swirl  causes  some 
streamline  spreading  and  central  pressure  reduction,  but  the  rotating  is 
too  weak  to  induce  recirculation.  For  S  >  0.6  or  so,  there  may  be  toroidal 
recirculation  zones,  but  the  flow  in  the  immediate  vicinity  of  the  axis  of 
rotation  may  be  downwind  at  all  axial  positions.  For  S  >  1  ,  an  "equivalent- 
bluff-body",  or  bubble-type,  recirculation  zone  often  lies  about  the  axis 
of  rotation.  This  large-swirl  case  is  the  one  of  interest  here.  In 
furnaces,  the  "bubble"  can  be  open-ended  in  that  the  recirculation  zone 
extends  from  near  the  inlet  down  the  entire  axial  length  of  the  chamber, 
with  the  result  being  a  nearly  columnar  flow  field  (Khalil,  El-Mahallawy 
and  Monei b  1 976) . 
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For  such  large-swirl  cases,  one  has  a  flow  that  is  well  described  in 
time  avarage  by  an  inviscid  incompressible  axisymmetric  swirling  model 
(Bossel  1973).  In  fact,  the  relevant  quasilinear  second-order  elliptical 
partial  differential  equation  is  precisely  that  which  is  there  believed 
pertinent  to  the  two-cell  severe  atmospheric  vortex.  The  vortex-breakdown 
analogy  between  severe-atmospheric-vortex  flows  and  recircul atory  combustor 
flows  seems  more  than  a  possible  superficial  resemblance  to  be  noted  in 
passing.  The  equation  is  usually  expressed  in  cylindrical  polar  coordinates, 
and  involves  both  axial  and  radial  partial  derivatives,  and  is  not  fully 
specified  until  flow  conditions  pertinent  at  the  inlet  and  at  the  outlet 
are  specified  (leibovich  1978). 

That  the  mathematical  formulation  of  the  strong-swirl -stabil  ized 
furnace  flow  is  so  similar  to  that  of  the  two-cell  atmospheric  vortex 
is,  in  restrospect,  obvious,  because  the  physical  phenomena  have  appreciable 
similarities  (Carrier,  Fendell  and  Feldman  1980).  In  both  instances,  there 
is  an  annular  flow  displaced  from  the  axis  of  symmetry,  in  which  rapidly 
swirling  flow  moves  axially  "downwind"  (the  "eyewall").  Enveloping  the 
axis  there  is  a  "filling"  flow  whose  axial  direction  of  movement  is  partly 
reversed  to  that  of  the  annular  fluid.  This  reversed-flow  fluid  in  the 
"bubble"  (or  "eye")  may  be  swirling,  but  the  principal  interest  is  in  (1) 
ascertaining  the  location  of  the  interface  distinguishing  annular  (eyewall) 
fluid  from  bubble  (eye)  fluid,  and  (2)  details  of  the  annular  flow  outside 
the  bubble.  (What  flow  there  is  inside  the  bubble  may  be  intricate  and  it 
is  probably  not  critical  to  ascertain  the. details  for  many  purposes.) 


It  is  not  the  objective  here  to  deal  exhaustively  with  similarities 
(and  also  distinctions)  between  recirculatory  bubble-type  swirling  furnace 
flows  and  two-cell  severe  atmospheric  vortical  storms.  Nevertheless, 
several  points  may  be  added.  First,  the  outer  boundary  on  the  annular 
flow  is  prescribed  in  the  case  of  a  furnace:  in  general  it  is  the  impervious 
outer  solid  wall  of  the  ini et-diffuser-dump  geometry.  For  the  atmospheric 
vortex,  there  is  no  lateral  surface,  and  a  second  free  bounding  surface 
must  be  identified  in  the  course  of.  solution:  it  is  the  (nonsolid)  boundary 
between  the  eyewall  and  outer  vortex  of  the  hurricane.  Thus,  the  atmospheric- 
vortex  case  has  a  free  surface  bounding  the  eyewall  on  both  its  near-axis 
and  outer  limits;  the  atmospheric-vortex  case  is  more  difficult  in  this 


regard.  Second,  whereas  the  swirling  in  the  annulus  is  often  of  rigid-body 
type  in  a  combustor,  an  analytically  simple  special  form  that  reduces  the 
quasilinear  partial  differential  equation  to  a  far  more  tractable  linear 
equation,  the  swirling  in  the  annulus  in  the  atmospheric  vortex  is  likely  to 
be  of  a  more  complicated,  less  tractable  nature.  Nevertheless,  the  total 
pressure  head  is  invariant  from  one  streamline  to  another  in  the  atmospheric 
vortex,  while  it  is  almost  certainly  not  invariant  in  the  furnace.  While 
these  characteristics  are  somewhat  off-setting  in  the  complication  that  they 
entail,  the  atmospheric  vortex  case  ends  up  the  more  challenging.  Third, 
in  both  flows  there  are  heat  releases:  from  chemical  exothermic i ty  in  the 
furnace,  and  from  the  latent  heat  of  condensation  in  the  hurricane.  In 
neither  case  is  this  heat  release  by  itself  important  enough  to  necessitate 
introduction  into  the  dynamics  of  compressibility  for  a  physically  useful 
quantif ication .  However,  the  atmospheric  vortex  extends  from  ground  level 
to  tropopause,  and  over  this  scale  of  altitude  compressibility  enters; 
of  course,  a  furnace  certainly  has  no  comparable  extent.  Nevertheless,  if 
one  confines  attention  to  just  the  very-near-ground  portion  of  the  eye/eyewall 
interface  in  a  two-cell  vortex  (a  portion  referred  to  as  the  turnaround  or 
corner,  in  that  low-level  swirling  influx  becomes  swirling  upflux),  incom¬ 
pressibility  suffices  for  the  atmospheric  vortex  as  well.  Fourth,  one  of 
the  challenging  matters  in  the  elliptic  flows  that  characterize  bubble-type 
furnace  flows  and  atmospheric  vortex  flows  is  appropriate  characterization 
of  the  "outlet"  flow.  Elliptic  flows  require  specification  of  flow  constraints 
at  all  "boundaries";  such  specifications  at  the  inlet  and  the  lateral 
"boundaries"  are  straightforward  to  state  physically  even  if  sometimes 
awkward  to  manipulate  mathematically.  However,  one  has  less  experience 
concerning  what  is  appropriately  specified  at  the  exit  plane,  such  that  the 
formulation  is  neither  underspecified  or  overconstrained .  It  is  to  be 
remembered  that  the  exit  plane  of  the  annular  fluid  may  be  the  entrance 
plane  of  the  core  ("bubble"  or  "eye")  fluid,  and  the  mathematics  must  be 
"informed"  about  the  properites  of  such  mass  as  enters  the  computational 
domain.  In  the  case  of  a  hurricane  it  is  known  that  the  well -developed 
eye  extends  from  tropopause  to  virtually  sea  level  and  contains  virtually 
nonrotating  air,  whereas  an  open-ended  reversed-flow  domain  in  a  swirling 
furnace  may  have  appreciable  swirl  within  the  air  entering  at  the  downwind 
extremity. 
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2.  A  BOUNDARY-VALUE  PROBLEM  FOR  THE  TURNAROUND 

An  inviscid  incompressible  steady  axisymmetric  model  of  the  turn¬ 
around  region  in  a  severe  vertical  vortex  is  adopted,  in  cylindrical 
polar  coordinates  with  origin  at  the  ground  on  the  axis  of  rotation  (and 
of  symmetry)  (Fig.  1 ) • 

The  azimuthal,  radial,  and  axial  velocity  components  (v*,  u*.  and  w*, 
respectively),  in  view  of  conservation  of  angular  momentum  and  of  conser¬ 
vation  of  mass,  may  be  written  [r*  =  r£V*  ,  where  r*,  is  the  (given)  radius 
of  peak  swirl  in  region  I,  and  where  V*  is  the  (given)  value  of  that  peak 
swirl] 


(r*v*)2  =  r*2F(ii<); 

(1) 

_  r*h*  9<|/  . 

•r  u  -  *  az*  ; 

(2) 

r*w*  =  r*h* 
r  W  tt  3r*  • 

(3) 

Super  asterisk  denotes  a  dimensional  quantity;  no  asterisk,  a  dimensionless 
quantity.  The  cylindrical  radial  coordinate  is  r*;  the  axial  coordinate,  z*; 
the  maximum  height  of  the  surface  inflow  layer,  h*;  the  streamfunction  for 
the  secondary  flow  (involving  velocity  components  u*,  w*.  only), the 
function  giving  distribution  of  angular  momentum  with  streamfunction,  F. 


Conservation  of  radial  momentum  is  (p£  is  density) 


u* 


w* 


3U*  _3_ 

dz*  dr * 


r*2FU) 
„*3  ' 


(4) 


Bernoulli's  equation  is 


*  g*z* 

2r*^ 


(5) 
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in  (5)  the  Bernoulli  constant  is  taken  as  universal  for  all  streamlines 
emerging  from  region  II  in  view  of  the  known  relation 


[u*(r*,z*)]2 
- 2 - 


+ 


(6) 


and  the  fact  that  the  pressure  above  hydrostatic  is  axially  invariant  for 
z*<h*,  within  the  conventional  boundary-layer  approximation.  Within 
the  cyclostrophic  approximation. 


P*  -  pf  =  (PJ/2 }V*2  =  (p^/2) (r*2/r^) ,  (7) 

where  p*  is  the  ground-level  ambient  pressure  (given). 

From  a  crude  treatment  of  the  turnaround,  to  be  published  elsewhere, 

(h*/r^)  =  0.2  for  relevant  values  of  r^,  V*.  where  typical  (given)  values,  for  a 
severe  tornado,  are  r*  =  160  m,  V*  =  100  m/s.  In  any  case,  h*  and  r*  are 
taken  as  known.  Thus,  in  addition  to  its  use  in  (2)  and  (3),  h*  is  employed 
to  nondimensionalize  the  independent  variables: 

x  =  Trr*/h*  ,  y  =  iTZ*/h*.  (8) 

Subtraction  of  the  radial  derivative  of  (5)  from  (4)  gives 

*  (l  ♦,),  *  *yy  *  ^  ■  0  •  (9> 

where  subscripts  x  and  y  denote  partial  differentiation.  The  swirl  profile 
F(ij>)  is  known  from  conditions  holding  at  x  *  x1  =  irrf/h*  =  rr/ h.  It  is  known 
from  published  solution  of  region  II  that  F  is  a  monotonically  increasing  function 
of  ip,  such  that  F(0)  =  0  and  F(I)  =  1,  where  the  datum  ( x >x ^ , 0 )  =0  is  adopted. 
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and  from  (2)  and  (8), 


’Kxj.n)  = 


(10) 


where  u  3  u*/V*,  and,  for  completeness,  v  3  v*/V*.  Clearly,  for  F  linear 
in  (9)  becomes  Poisson's  equation,  and  for  F  quadratic  in  <h,  (9) 
becomes  Helmholtz's  equation.  Although,  from  known  results  for  region  II, 
such  simple  forms  for  F(ii>)  are  not  detailed  replications  for  (r*v*)/T*  of 
interest,  still  their  tractability  urges  their  adoption. 

The  pressure  along  the  streamline  that  separates  from  y  =  0  at  x  = 
remains  at  p^,  for  the  model  of  a  nonrotating  "eye"  isobaric  at  altitudes 
at  which  gravity  plays  no  role.  Then,  from  (5),  since  F(0)  =  0, 

2 

(*/  +  Uy)2  =  on  ^  =  0  .  (II) 

The  pressure  along  the  streamsurface  =  I,  which  passes  through  the 
circle  x  =x1,y  =  it,  should  be  consistent  with  the  pressure  in  the  potential- 
vortex  region;  hence 


P(x,y)  =  Pj  + 


(12) 


from  the  cyclostrophic  balance,  where  the  nondimensional ization 


P 


P* 


1  n*\J*2 

2  dv 


(13) 


has  been  adopted.  From  (5),  (12),  and  (13),  with  neglect  of  the  gravitational 
term,  since  F(I)  =  1, 


<h 


2 

x 


0  on  i//  *  I  . 


(14) 
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This  condition  merely  states  that  there  is  no  contribution  from  the  secondary- 
flow  velocity  components  u*.  w*  to  the  pressure  field,  where  region  III 
interfaces  with  region  I. 

Equations  (9),  (11),  and  (14)  constitute  the  boundary-value  problem 
of  interest,  when  supplemented  by  a  relation  for  F(i|0,  obtained  from  pre¬ 
viously  executed  analysis  of  region  II.  By  (1),  (2),  (6),  and  the  fact 
that  w(x^,y)  *  0,  the  relation  for  F(ip)  gives  implicitly  an  expression  for 
Y(x-|,y),  0  <  y  <  it,  with  which  to  initiate  the  analysis  (see  below).  Finally, 
it  is  anticipated  that  ip(x ,y ) ,  as  given  by  the  boundary-value  problem  is 
periodic  in  y;  this  periodicity  reflects  the  fact  that  the  formulation  must 
be  revised  after  completion  of  about  one  cycle.  More  explicitly,  it  is 
anticipated  that  the  locus  of  any  streamline,  say  ip  =  1/2,  as  y  increases, 
decreases  to  a  minimum  value  of  x,  then  increases  in  x  to  recover  its 
initial  value  x  =  x^  and  to  achieve  a  peak  value  in  x,  before  decreasing 
back  to  x  =  x1  to  begin  another  period.  It  is  to  be  expected  that  the 
undulating  flow  so  generated  will  be  unstable  somewhere  along  its  trajectory, 
but  our  present  job  is  to  find  that  flow  and  the  stability  problem  is 
deferred  until  that  is  done.  Of  particular  interest  are  the  amount  by 
which  v  exceeds  the  value  associated  with  x  5  Xj,  (/  8  I,  and  the  values  of 
x  and  \p  at  which  the  maximum  occurs. 

For  explicitness,  a  rough  characterization  of  the  function  iHx-j.y), 

0  <  y  <  7i,  is  now  presented.  In  the  largely  inviscid  portion  of  region  II, 

(6)  may  be  rewritten  as 

u2(x,y)  +  v2(x,y)  =  (x^x)2  .  (15) 

From  (1),  (2),  and  (15),  one  obtains 

-u(xry) - ji —  =  {l  -  [v(x1,y)]2J  =  [1  -  F(*)]1/2,  (16) 
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where  C  <_  F ( ^ )  <_  1  for  0  <_  ^<_I.  Thus,  'f'(x^.y)  is  given  by 

♦(x, ,y) 

0 

where  i|»  I  as  y  -*•  *  (cf.  {10)).  For  F(i^)  =  1  * 

*  *  y(l-  &)  *  1  *  \  *  y2(xi*y) =  ¥  (*  •  h)  •  (18) 

For  F(*)  =  u/l)2, 

*  =  21/2sin(y/2),  I  =  21/2*> v^Xj.y)  *  sin2  (y/2)  .  (19) 

In  fact, 

v(x1,y)  *  1  -  exp  (-  j  y)  ,  c  =  3  ,  (20) 

might  be  more  realistic,  but  is  far  less  tractable. 

In  summary,  if 

*x,y)  =  .  G'fT)  =  ^  .  (21) 

Then  the  boundary-value  problem  is,  for  0  f.  y  ^  1, 

X(^4'x)  +  fyy  +  G ’{')')  =  o  ;  (22) 

f 2  +  f2  on  f  =  0,  'fx  +fy  a  0  on  t  *  1;  (23a) 
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(23b) 


,y)  given,  with  o  periodic  in  y. 

One  might  consider  interchanging  the  roles  of  the  dependent 
variable  i>  and  the  independent  variable  x,  and  using  perhaps  some  technique 
within  the  framework  of  the  method  of  weighted  residuals,  perhaps  with 
multigrid  procedures  (Brandt  1977),  to  extract  the  desired  information 
from  (21)  -  (23).  However,  first,  alternative  procedures  are  used  to  study 
(21)  -  (23). 

3.  TREATMENT  OF  A  SEPARATED  LAYER  WITHOUT  STRUCTURE 


A  momentum  balance  in  an  axisymmetric  separated  boundary  layer  without 
structure  is  examined.  The  thin  sheet  (or  "eye  wall")  is  the  demarcation 
between  any  "eye,"  isobaric  at  pressure  p^,  and  a  potential  vortex  with 
radially  dependent  pressure  given  by  (12).  At  r*  *  r*.  p*  in  the  sheet  equals 
that  in  the  "eye"  (Fig.  2). 

Henceforth  in  this  section  the  symbol  r*(z*)  denotes  the  outside 
surface  of  the  "eye  wall"  without  structure.  At  the  point  A,  given  by 
[r*(z*)»  z*]  ,  the  principal  radii  of  curvature  are  |l  +  [r*' (z*)]^|^/r*"(z*)  , 
the  radius  of  curvature  in  a  plane  containing  the  axis  of  symmetry  and  the  stream¬ 
line  on  which  point  A  lies,  and  r*(z*)|l  +[r*'  (z*)]^|^,  the  radius  of 
curvature  in  a  plane  perpendicular  to  the  streamline  on  which  point  A  lies. 

Super  prime  denotes  ordinary  derivative  with  respect  to  the  argument  of  the 
function.  The  velocity  component  in  a  plane  containing  the  axis  of  symmetry 
is  denoted  q*,  while  the  velocity  component  (swirl)  in  a  plane  perpendicular 

A 

to  the  axis  of  symmetry  is  denoted  v*.  Thus,  if  t  is  a  unit  vector  in  the 

A  AV 

plane  containing  the  axis,  and  0  is  perpendicular  to  t  and  refers  to  the 
azimuthal  component  in  a  cyl indrical-polar-coordinate  system,  then 

v*  =  q*t  +  v*0.  (24) 
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The  force  balance  perpendicular  to  the  thin  sheet  at  [r*(z*),  z*] 
equates  the  pressure  gradient  consistent  with  a  potential  vortex  outside 
the  sheet,  to  the  component  of  acceleration  perpendicular  to  the  sheet. 
There  are  two  contributions  to  the  acceleration,  each  involving  the 
square  of  a  velocity  component  over  an  appropriate  radius  of  curvature: 

where  h*  represents  a  coordinate  running  across  the  sheet. 

Variations  in  velocity  occurring  across  the  sheet  are  not  resolved 
(F*  ^  r*V*): 


2 

q*  r*" 


[1  +  (r*')2]3/Z  r*[l  +  (r*‘ )  ] 


r*W 


The  potential- vortex  form  of  v*  is  used  in  the  definition  of  A*. 


From  these  definitions  and  from  (12),  to  an  accuracy  that  serves 
current  purposes. 


A. ill),  _ 

'  r*  '  r*3fl  f  (r* 1 


3*y*2r*» 


[l  *  (r*1)2]  [l  +  (r*‘  )2J' 


The  following  nondimensionalization  is  introduced: 


r*  —  _  z*  _  2A*  _  2B* 

rl  •  rl  rl  rl 


If  A*  =  (h*/2),  B*  =  (h*/2),  seemingly  reasonable  values,  then 


„  i  h*  _  . 
a  -  pr  =  n  , 

rl 


„  .  h* 

8  =  pr  =  h 

rl 


Under  (29),  (28)  becomes  (a, 6  specified) 


i?  *  O  ?  1/2  "  _  ;  3/2 

x  XJ(1  +  x,z)  (1  +  x,Z) 


_  a  1  +  x’2  1  /x2  -  l\  ,  71 2 , 3/2 

(  -  S  B  \  ^2  /  (1  +  x  >  • 


In  this  translationally  Invariant  equation,  it  is  taken  that  7(7  »  0)  =  1. 
It  Is  shown  below  that  7  is  periodic  in  7.  Solution  is  sought  for  positive 
and  negative  z,  where  the  boundary  conditions  are 


x (0)  *  1 


(32) 


x'(0)  *  x^,  given  const.  >  0. 


(33) 


The  boundary  conditions  preclude  odd  or  even  solution  for  7  in  z.  Sought 
are  7  +  ,  the  largest  value  of  7,  which  occurs  where  7*  *  0,  and  7'  ,  the 
smallest  value  of  7,  which  also  occurs  where  7‘  =  0.  In  that  7~  <  l, 
there  is  "overshoot"  of  the  swirl,  and,  from  (12),  there  is  decrease  of 
pressure  from  the  value  that  holds  in  the  "eye",  i.e.,  in  0  <  r*  <  r*(z*); 
it  is  reiterated  that  the  magnitude  of  the  swirl  overshoot  and  pressure 
decrease  are  of  particular  interest. 

While  numerical  integration  of  (31)  -  (33)  is  required  ultimately, 
some  preliminary  treatment  is  helpful.  Since  (31)  is  translationally 
invariant  in  z,  phase-plane  analysis  is  introduced  (Fig.  3  ); 

*  F(x)  *>  =  — ffl  -  (34 

dz  d !r  dz  dz  dx  dx 

Thus , 


bpp1 


377 


(1  +  P2) 


-3  -2  1/2 

xJ  (i  +  r ) 


(35) 


The  slope  is  infinite  at  F  =  0  and  is  zero  where 


j»a 


+  P2) 


1/2 


1 

? 


-1  =  0 


(36) 


the  intersection  of  the  curve  of  (36)  with  the  P  =  0  ray  is  given  by 


=ar +  ^ "  i =  0  * 


(37) 


where  x*  =  1  for  a  =  0,  x*  >  1  for  a  >  0.  For  1  >>  a  >  0, 
x*=  1+  (a/2)  -(3a2/8)  +  .... 

Although  the  following  development  is  not  pursued  to  the  extent  of 
obtaining  results,  it  may  be  worth  noting  that  if 


then  (35)  becomes 


H(x)  =  (1  +  P2)'172  , 


(38) 


-6H' 


aH  .  i  , 
=3  +  Z?  -  1 

X  X 


(39) 


If 


T  =  X  '2  , 


(40) 


then 


2b£-  aH  =  x-l/2  .  T-3/2 


(41) 


or 


•  (zk)I/!  i1  *  f) (if)  erfc  (s)'/2 


+  E  exp 


(w)’ 


(42) 
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where  E  is  a  const,  of  integration.  From  (38),  (40),  and  the  definition 
"F  =  (dx/dz),  one  may  write  formally 

—  =  {1  +  [h(t)]'211/2  =  d<T ~l/-l  . 

dz  '  ’  dz 

Thus, 


dz 


{i  *  CHh)r*}w ' 


(44) 


where  H(t)  is  given  by  (42). 

For  the  special  case  a  =  0,  i.e.,  no  swirl,  multiplication  of  (35) 
by  x'  yields 


s .  .  Lx  -  ir  +  s 

(1*  x  "•  (45) 

where  the  constant  of  integration  has  been  written  as  (6/m),  with 


m  =  (1  +  x^2)1/2  (46) 

for  consistency, 
finite  constant, 
so  (1  +  7;2)1/2 
At  x'  *  0, 

5-4  .  1  +  |  (1  -  *  [fd  -  i)]1/Z  [z  +  !  (1  -  i)]l/Z  (47) 

For  x^-^0  so  m  1,  x  ±merge  to  unity;  this  case  involves  a  vertically 

separating  surface  inflow  layer,  and  hence  no  overshoot.  For  7' -►  ®  so 

±  0 

m  ■*  «,  x  remain  bounded: 


It  is  recalled  that  xQ'(  =  dx(0)/dz)  is  a  given  positive 
Inspection  of  (45)  reveals  that  x'  is  maximum  at  7*1, 
is  the  maximum  value  of  (1  +  x'  )  ,  whence  the  symbol  m. 
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this  case  involves  effectively  horizontal  inflow  of  the  separating  surface 

inflow  layer,  and  leads  to  the  minimum  value  of  x  "  for  fixed  3.  For 

3  =  0.1,  x+  =  1.37  and  x'  =  0.73;  for  3  =  0.2,  x+  =  1.56  and  x‘  =  0.64. 

The  plausible  range  of  6  =  0(0.2),  from  a  crude  analysis  of  the  turnaround 

to  be  published  separately;  it  is  recalled  that  £  is  associated  with  in-plane 

motion.  These  results  suggest  what  proves  to  be  a  general  trend:  x" 

decreases  as  x'  and  8  increases, 
o 

Finite  values  for  a  indicate  finite  swirl;  increasing  ;  yields  larger  x* 
and  smaller  overshoot.  From  results  of  numerical  integration,  for  finite  j., 
x"  approach  finite  values  as  P  •-  °°  for  fixed  S.  For  the  plausible  values 
a  =  3  =  0.2,  for  x^  =  2,  5,  10,  the  corresponding  values  of  7~  =  0.77, 

0.72  ,  0.70.  Hence,  swirl  overshoots  in  the  range  of  about  10°I  seem  plausible, 
but  no  more;  estimates  that  swirl  in  the  turnaround  exceeds  the  swirl  at 
x  =  1  ,  z  =  h  (in  terms  of  the  definitions  of  (29))  by  about  lOO'I  (Lewellen  1977 
are  excessive  according  to  this  analysis.  Further  results  are  given  in 
Table  1  and  graphical  presentation ' i s  given  in  Figures  4  and  5. 

4.  REFORMATION  OF  THE  BOUNDARY-VALUE  PROBLEM 

A  change  of  variables  from  (r*,z*)  to  (^*,G),  where  the  angle  i  gives 
the  orientation  in  the  (r*,z*)  plane  of  the  streamline  relative  to  a  line 
parallel  to  z*  =  0  (Fig.  6),  is  developed  for  the  boundary-value  problem 
introduced  in  Section  2.  The  restriction  to  constant-density  flow  is 
rel axed . 

The  streamfunction  '+*  is  redefined  by 

,  .  .  .  _ 

P  q  r*  sine  =  ^  , 

p*q*r*  cos©  =  . 

Here  r*  again  denotes  the  conventional  cylindrical  radial  coordinate  at 
which  the  velocity  is  to  be  calculated,  q*  is  the  velocity  in  the  (r*,z*) 
plane,  and  n*  (introduced  below)  is  the  distance  normal  to  the  streamline. 


(49) 

(50) 
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The  quantity  q*  sinO  is  the  vertical  velocity  w*,  and  the  quantity  -q*  cos 
is  the  radial  velocity  u*. 

The  change  of  variables  from  (r*,z*)  to  (v*,3)  requires  that  both 
(3/3r*)  *  and  (s/sz*)  *  be  expressed  in  terms  of  (3/3y,*)  and  (3/33)  *, 
where  the  subscript  is  held  fixed.  Consider  a  function  f*  that  depends 
on  r*  and  z*,  and,  therefore,  on  ^*  and  0.  One  may  write 


df*  = 


■  m, * 

■  (»),  [M„  *•  •  (»)  «■]  *  - 

=  P*q*r*  sinG  dr*  +  p*q*r*  cos©  dz*  +  ^ 


Henceforth,  the  quantity  being  held  fixed  is  no  longer  explicitly  written, 
for  brevity. 


The  usual  transformation  rules  apply,  i.e., 


0  =  _ Q 

.1*  *  U 


O  :  *  _  O' 

J*  ’  v  r*  ’  J*  ’ 


where 


J*  =  r%*z*0  -  rV*z*  ‘ 


(52a  ,b,c  ,d) 
(52e) 


and  subscript  here  denotes  partial  differentiation. 

Division  of  (51)  by  dz*  and  letting  dz*  0  as  dr*  is  held  fixed  at  zero, 
one  obtains 


si mi larly , 


f*z*  =  f%*  p*c>*r*  cos0  +  f*Q  0z*  * 


f*  *  =  f*  p*q*r*  si n0  +  f*  0  * 
r*  ^*  ^  0  r* 


23 


From  (52),  (53),  and  (54), 


f*  *  =  p*q*r*  cosG  f*  *  +  -5 — -* - —  f*r  ; 

2  'v  r  p*z  0  '  r  Gz  v*  ° 


fV  -  p*p*r*  si"0  f>  ~  r-^V^rT-zTT  f*G  ’  (56) 

\p*  Q  o  p* 

In  particular,  the  directional  derivative  j_n  the  r*  ,z*  frame,  with  direction 
normal  to  a  constant  ij;*  line,  is 


*n*  -  f*r*  sinO  +  f*2*  cosG  ; 


i  .e. , 


r*i(*  cosG  -  z*.  *  sinQ 

f*  =  n*n*r*  f*  +  _ id _ , _ '£ _  f*  l 

f  n*  p  q  r  f  ip*  r*  *z*  -  r*rz*lir  r  0  * 

Y  ip*  0  G  p* 

If  one  adds  (55)  with  f*  replaced  by  z*  to  (56)  with  f*  replaced  by  r*, 
one  gets 

z*  cosG  +  r*  sinG  =  — rL -»  .  ( 

P*  p*  p*q*r’r 


The  quotient  (4d)  t  (4c)  gives 

z*0  =  "  r*0  tan0  • 

Also,  by  use  of  (58)  and  (59)  to  eliminate  z*  in  the  expression  for  J* 
given  in  (52e), 


J*  =  -r*,*  r*n  tan©  -  r*r,(  * 

'I'  0  S\P  q*r 


0lp*q*r*  cosG  ^ 


r*  *  tanG 


p*q*r*  cosG 


Furthermore, 


r%*  cos0  -  2>  sinG  =  rV  0050  ■  [p*q*r*  cos©  -  rV  tan0J  sin0 

r* 

.  p*  _  tanG 
cosG  "  p*q*r*  * 


Substitution  of  (60)  and  (61)  in  (57)  gives 


tanO 


f*n#  =  c*q*r*f*^  -  c?^.  ^  q  r  p*q*r*  cosQ 


=  D*q*r*f*  __  s i  nO 

4*  n  „  o  r  ~ 


0 


•j 


0 


Thus,  the  conservation  of  momentum  in  the  direction  normal  to  a  surface 
G*  =  const,  is  given  by 

r* 

p*qV*p*  .  p*q*r*  p*  +  p* 

V  r  Q  u  r  Q 

2  2 

p*v*  sinO  ,  o*q*  cosO  *  * 

=  t — -j - +  — ** - P  g*  cos-  . 

0 

[The  middle  term  on  the  right-hand  side  is  obtained  thus: 

0  q  3s*  D  q  fas*/30)  0  q  ( 3r*/50)  ’ 

where  s*  is  distance  along  a  streamline.  The  principal  curvatures  of  a 
surface  <p*  =  const,  in  the  (r*,z*)  plane  are  thus  seen  to  be  (3G/3S*)  and 
sinG/r*.]  Rearrangement  gives 


(pt‘ 


2 

v*  sine  .  g*cos0 

^  •  j.  j. 


*  2 
r*  q* 


*  cosO 
r* 


q*r* 


)r*0  '  (rV  -  p'q'r* )  j 


(62) 


If  one  divides  (58)  by  cosG,  differentiates  the  result  with  respect  to  G, 
and  then  subtracts  from  this  result  the  derivative  of  (59)  with  respect 
to  0*,  one  obtains 


r*  =  cos  01 


( _ I _ \ 

\p*q*r*  cos0/Q  ‘ 


(63) 
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The  pertinent  equations  include  an  equation  of  state,  a  (pseudo) 
adiabatic  relation  ,  conservation  of  angular  momentum 

r*v*  =  F*(ip* )  ,  (64) 

and  Bernoulli's  equation 

+  ^2“  +  -p  +  9*z*  =  G*(v*)  ,  (65) 

where  F*(tp*),  G*(^*)  are  implied  by  the  boundary  conditions  and  are  taken 
as  specified.  Equations  (59),  (62),  and  (63)  supplement  this  list.  There 
are  seven  equations  for  the  seven  unknowns  q*,  v*.  r*,  z* ,  p*,  p*.  and  T*. 
Each  unknown  is  a  function  of  < p*  and  0. 

The  appropriate  initial  conditions  are 

p*  •  ,  r*  =  r^  ,  F*(ip*)  given,  G*(ip*)  given  on  0  =  GQ  =  0, 

0  <  **  <  ^  .  (66) 
At  0  *  0  =  0,  0  <  ip*  <  ,  one  anticipates,  for  r*  =  r* 

0  maX  I 


2  2  •  2 
q*  +  v*£  =  V*  , 

(67a) 

pf  -  p*  =  o*g*z*  . 

The  boundary  conditions  for  small  z*  are 

(67b) 

p*( <P*  a  0,0)  =  p*  , 

(68) 

The  reference  here  is  to  a  relation  among  the  thermodynamic  variables 
obtained  from  the  second  law  of  thermodynamics .  In  the  absence  of  conden¬ 
sation  of  water  vapor,  the  dry-adiabatic  relation  suffices. 
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All  quantities  with  subscript  unity  are  given  constants,  C-*  is  calculable, 

iTIa  o 

and  r*  =  r^  V*,  with  V*  given. 

In  Figure  6,  p*  is  denoted  o*d*.  and  the  possibility  that  the  boundary- 
value  problem  possesses  a  reflection-type  solution  at  low  altitudes  is 
explored  schematically. 

With  the  boundary-value  problem  now  displayed  in  (o*,0)  coordinates, 
it  is  moot  whether  this  formulation  is  preferable  to  that  in  (r*,z*) 
coordinates. 

5.  GENERALIZATION  TO  A  VARIABLE-DENSITY  MODEL— FURTHER  COMMENTS 

For  a  formulation  valid  through  the  depth  of  the  tropopause,  it  is 
useful  to  review  tenets  of  the  modeling  at  this  point,  since  greater  detail 
concerning  the  thermodynamics  is  to  be  added. 

When  the  state  of  the  ambient  atmosphere  is  known,  and  when  the 
circulation  7*  of  the  potential  motion  is  given,  one  can  calculate  the 
inviscid,  quasisteady  motion  which  is  consistent  with  any  specified  central- 
core  pressure  distribution,  the  logical  extremes  of  which  are  (1)  the  dry- 
adiabatic  condition  consistent  with  the  total  enthalpy  of  the  ambient  air 
at  sea-level,  and  (2)  something  weaker  than  the  moist-adiabat  corresponding 
to  that  same  sea-level  condition. 

The  general  character  of  the  configuration  is  given  in  Figure  7.  In 
region  I  any  radial  and/or  vertical  motion  is  supposed  to  be  so  small  that 
it  is  ignorable  in  the  dynamic  balance.  The  angular  momentum  7*  is  pre¬ 
scribed,  most  simply  as  a  constant;  ultimately,  interest  centers  on  the 
case  where  7*  is  a  function  of  z*.  Region  II  is  the  boundary  layer  which 
transports  swirling  fluid  inward;  r*  is  the  radius  at  the  pressure  in  the 
potential  vortex  at  ground  vortex  is  the  same  as  the  prescribed  pressure 
at  r*  =  0,  z*  =  0.  Region  III  contains  updraft  fluid  and  its  state  trajec¬ 
tory  is  dynamically  moist  adiabatic;  i.e.,  the  total  enthalpy  is  constant. 
Region  IV  is  nearly  stagnant  and  contains  either  updraft  air  (nearly  moist 

adiabatic)  or  recompressed  air  from  z *,  the  altitude  at  which,  in  the 

0 

ambient,  the  total  enthalpy  is  the  same  as  that  at  ground  level.  Typically, 
that  ambient  enthalpy  is  of  the  character  depicted  in  Figure  8. 


It  has  already  been  noted  in  Section  4  that  the  boundaries  of  region 
III  are  not  known;  the  coordinates  r*  and  z*  are  not  particularly  well 
suited  to  a  successful  analysis.  On  the  other  hand,  if  one  defines  a 
streamfunction  ip*  and  a  "velocity  direction  0",  the  boundaries  are  fixed 
in  the  (-^*,0)  plane  and  it  may  be  useful  to  formulate  the  dynamic  balances 
in  that  framework. 

In  writing  the  equations  for  the  conservation  of  momentum  for  the 
inviscid  flow  of  an  air,  water-vapor  mixture,  one  usefully  adopts  the 
hypothesis  that  no  supercooling  occurs,  but  that  once  the  temperature 
drops  far  enough  to  imply  saturation,  water  condenses  at  precisely  the 
rate  needed  to  keep  the  vapor  at  the  equilibrium  saturation  level. 

The  angular  momentum  of  each  particle  is  conserved,  and  the  total 
head  is  preserved  on  each  streamline.  Momentum  conservation  normal  to  a 
streamline  and  geometric  relations  that  supplement  these  statements  were 
derived  in  Section  4. 

Finding  the  locus  of  thermodynamic  states  that  characterize  the  moist 

adiabat  based  on  sea-level  ambient  conditions;  assignment  of  an  altitude  to  the 

tropopause,  denoted  z*  in  Fig.  8;  and  determination  of  the  pressure  distri- 

o 

bution  with  height  in  region  IV,  denoted  PQ(z*)  in  Fig.  7  --  these  steps 
are  standard  (Fendell  1974). 

Finding  the  pressure  in  region  I  compatible  with  the  ambient  atmos¬ 
phere,  p*(r*,z*),  is  now  discussed.  If  one  discards  the  subscript  unity. 


P*r*  =  P*v*2/r*  , 


P*z*  =  -pV 


hence , 


v* 

1 _  D* 

r*  P  z* 


r* 


If  one  seeks  solution  in  the  form 


p*(n*) 


3*  f 


p*  m*(z*) 
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one  obtains  (if  prime  denotes  the  ordinary  derivative) 


P*tn*^g*s*'(r*)r*3  -  r*2(z*)m*'(z*)J  =  0  ; 


i.e.,  without  loss  of  generality, 

s *(r*)  = 


1 


2g*r* 

.*(!*)  ■  r  -r—  • 

■>0  r'Mz' ) 

For  the  special  case  of  a  vortex  in  region  I  that  is  invariant  with 
altitude,  F*  =  const. , 


m*(z*)  =  z*/r*2  ; 


the  choice  of  p*(n*)  which  matches  the  ambient  atmosphere  at  r* 


<< 


is 


P*  = 


pambient 


+  r* 


[In  this  approximation  the  ambient  formally  is  taken  to  hold  at  an  infinite 

radial  distance,  but  the  discrepancy  from  a  match  to  an  ambient  at  a  value 

of  (r*/r*)  of  only  ten  results  in  an  error  of  about  merely  one  percent.] 
o  l 

For  the  somewhat  more  plausible  case  of  a  vortex,  the  angular  momentum 
of  which  decreases  linearly  from  a  value  of  r*  at  sea-level  to  zero  at  the 
tropopause  z*  =  z*  , 


n*  = 


z*  - 
0 


r*Z 


2g*r* 


such  that  at  r*  -*■  ®, 


z* 


z*n* 

_ o  . 

n*  +  l  ' 


Thus  a  solution  of  the  first-order  partial  differential  equation  p*(n*), 
any  function  of  n*.  also  satisfies  the  boundary  condition. 
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P*(n*) — >  Pambient^2*^  ^or  r*  — >  one  takes  as  the  function 

~*2 , 


p*(n*)  ■*  p*  .  .  * 

K  v  ^ambient 


z*  + 


o‘(zo  ‘  z*} 


2q*r*22n 


1  + 


r0*2u;  -  2*) 
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*r *  z* 


In  region  III  the  pressure  on  the  streamline  contiguous  to 
region  IV  must  match  the  specified  p*(z*)  in  the  core,  and  the  pressure 
on  the  streamline  contiguous  to  region  I  must  match  the  pressure  in  the 
potential  vortex  Pgmbient^m*  -  s * ( r* ) ) ] -  These  statements  consti¬ 

tute  the  boundary  conditions.  At  r  =  r*,  0  <  z*  <  h*.  the  initial  value 
of  each  dependent  variable  is  specified  on  the  basis  of  the  known  solution 
for  region  II. 

Obtaining  the  solution  of  this  formulation  of  "eye  wall"  with  struc¬ 
ture,  throughout  the  depth  of  the  troposphere,  is  the  next  goal  of  this 
on-going  investigation. 


Table  1.  Numerical  Results  for  Extrema  of  the  Turnaround,  from  Integration 
of  the  Initial-Value  Problem  (31)-(33). 


a 

a 

a/ a 

«*(o>  * 

(*  >,.o 

<r>c.lc 

0.1 

0.1 

0.0 

2.0 

0.791 

0.789 

1.264 

1.266 

0.0 

0.1 

0.0 

5.0 

0.754 

0.751 

1.326 

1.329 

0.0 

0.1 

0.0 

10.0 

0.742 

0.738 

1.348 

1.352 

0.0 

0.1 

0.1 

100.0 

0.732 

0.727 

1.366 

1.372 

0.0 

0.2 

0.0 

2.0 

0. 718 

0.714 

1.392 

1.397 

0.0 

0.2 

0.0 

5.0 

0.671 

0.664 

1  489 

1.497 

0.0 

0.2 

0.0 

10.0 

0.656 

0.647 

1.524 

1.533 

0.0 

0.2 

0.0 

100.0 

0.644 

0.632 

1.553 

1.565 

0.1 

0.1 

1.0 

2.0 

0.820 

0.789 

1.295 

1.266 

0.1 

0.1 

1.0 

5.0 

0.777 

0.751 

1.348 

1.329 

0.1 

0.1 

1.0 

10.0 

0.762 

0.738 

1.366 

1.352 

0.1 

0.1 

1.0 

100.0 

0.750 

0.727 

1.381 

1.372 

0.1 
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Schematic,  not  to  scale,  of  an  inviscid  model  of  the  turnaround 
region  of  a  severe  vertical  quasisteady  axisymmetric  vortex. 

The  "eye",  region  IV,  is  isobaric  at  pressure  pf,  over  the  vertical 
extent  of  interest  here,  where  p?  is  also  the  pressure  at  r**r?, 
z* » h*  (because  the  pressure  is  approximately  invariant  across 
the  surface  Inflow  layer  II).  The  "eye  wall",  region  III,  is  demar 
cated  by  two  streamsurfaces ,  if»*(r*,z*)  *  const,  (the  position  of 
each  to  be  determined),  encompassing  the  mass  efflux  from  region  II 
It  is  recalled  that  a  potential  vortex  holds  in  region  I. 


Figure  2.  Schematic  of  the  location  of  an  inviscid  "eye  wall  without 

structure"  demarcating  the  interface  between  an  isobaric  non¬ 
swirling  "eye"  at  pressure  p?  and  a  potential  vortex.  The  sheet 
representing  the  "eye  wall"  lias  thickness  h*  0;  in  the  turn¬ 
around  region.  Its  displacement  from  the  axis,  as  a  function  of 
height  above  the  ground  plane,  is  denoted  r*(z*). 
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Figure  3. 


Phase-plane  properties  of  equation  (35),  where  P  =  (dx/dz),  with 
x  the  dependent  variable  (x  >  0)  and  z  the  independent  variable. 
Isoclines  of  zero  and  infinite  slope  are  noted.  At  x  =  1,  the 
distance  from  the  axis  of  synmetry  at  which  the  surface  inflow  layer 
separates,  a^  finite  positive  slope  is  adopted.  For  a  >  0,  x*  >  1; 
for  a  *  0,  x*  *  1.  The  sketched  trajectory  (solution  curve)  is 
a  limit  cycle  (closed  curve  indicative  of  periodic  behavior).  The 
periodicity  is  not  of  physical  interest. 


34 


Figure  6.  This  schematic  of  the  inward,  then  outward,  excursion  of  the 

separated  inflow  layer  with  structure,  indicates  the  bounding  stream 
surfaces  ip*  =  0  and  4)*  =  C.  The  in-plane  radial-axial  flow 
speed  is  denoted  q*;  the  swirl  about  the  vertical  axis,  v*.  The 
inward  excursion  need  be  solved  for  streamfunction-incl ination  angle 
0  *  0  to  9  =  (ir/ 2),  since  the  solution  fore  =  (ir/2)  to  0  =  tt 
may  then  be  obtained  by  reflection.  An  analogous  statement  holds 
for  the  outward  excursion.  The  "eye"  is  isobaric  at  pressure  pJ, 
but  the  potential -vortex  pressure  varies  with  radial  position 
r*  as  noted. 


Figure  7.  Another  schematic  of  the  postulated  four-part  model  of  the  struc¬ 
ture  of  a  mature  intense  hurricane,  of  radial  extent  r£  and  axial 
extent  z£  (the  tropopause).  The  potential  vortex  in  region  I  is 
here  ascribed  modest  axial  variation,  such  that  the  associated 
pressure  field  is  denoted  p-|(r*,z*).  The  low-level  swirling 
influx  (region  II)  grows  to  thickness  h*,  and  separates  to  form 
the  "eye  wall"  (region  III);  the  flow  speed  independent  of 
swirling  is  denoted  q*,  and  the  inclination  of  the  secondary  (i.e 
radial  and  axial)  flow  relative  to  a  horizontal  plane  is  charac¬ 
terized  by  angle  0.  The  pressure  variation  in  the  "eye",  region 
I V ,  is  entirely  axial  and  is  denoted  p£(z*);  for  the  idealization 
of  a  completely  dry,  cloud-free,  fully  developed  eye,  this  varia¬ 
tion  is  determined  by  dry-adiabatic  compression  of  tropopause- 
level  air. 


p*(mb) 
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A  typical  profile  of  the  total  static  temperature  (H*/Cn)  over  the 
ocean  in  the  equatorial  trough  and  in  the  subtropics  (1200n  mi  from 
the  equatorial  trough),  as  a  function  of  pressure  p*  and  altitude 
above  sea-level  z*.  The  loci  of  thermodynamic  states  for  a  typical 
cumulus  (cu)  and  also  for  the  undiluted  core  of  a  cumulonimbus  (cb) 
are  also  noted;  the  locus  (cb)  is  that  appropriate  for  moist-adia¬ 
batic  ascent  of  sea-level  air.  The  dashed  curve  is  the  typical 
ambient  referred  to  in  the  text:  a  well-mixed  lowest  layer,  a 
midtropospheric  minimum,  and  a  recovery  of  sea-level  value  at  alti¬ 
tude  z*  (defined  here  to  be  the  tropopause,  and  typically  about 
50,000  ft).  The  quantity  H*  is  the  sum  of  the  static  enthalpy,  the 
latent  heat  equivalent  of  the  water  vapor  present,  and  the  gravita¬ 
tional  potential  energy;  the  quantity  c£  is  the  specific  heat  at 
constant  pressure,  which  may  be  taken  to  be  effectively  that  of  dry 
air. 
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